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Scattering cold particles on an TV-slit grating is shown to reproduce an interference pattern, that 
manifests itself in the near- field region as the fractal Talbot carpet. In the far-field region the 
pattern is transformed to an ordinary diffraction, where principal beams are partitioned from each 
other by (N — 2) weak ones. A probability density plot of the wave function, to be represented 
by a gaussian wavepacket, is calculated both in the near- field region and in the far-field one. 
Bohmian (geodesic) trajectories, to be calculated by a guidance equation, are superimposed on the 
probability density plot well enough. It means, that a particle, moving from a source to a detector, 
passes across the grating along a single bohmian trajectory through-passing one and only one slit. 
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I. INTRODUCTION. 

Wave interference is a most impressive phe- 
nomenon be it induced by waves on water, acous- 
tic waves, or electromagnetic ones - radio-waves, 
light waves, 7-radiation. Quantum-mechanical 
experiments dealing with interference phenom- 
ena p] display interesting interference phenom- 
ena. They involve self-affine fractal quantum 
evolution of the probability densities and quan- 
tum revivals [2J [3J [5J |6]. There is a special 
interest to the phenomena from the side of quan- 
tum computation and communication [7J El E] ■ 

Most fantastic interference effects are disclosed 
in near-field region, i.e., in the vicinity of an in- 
terference grating. In this region complex wave 
interference shows very exotic patterns named 
in literature as the Talbot carpets [21 |3j Hj. 
They disclose fractal-like self-similar structures. 
Henry Fox Talbot was first who observed in 1836 
such an effect in a near-field region [32]. Fig. [TJ 
for example, demonstrates the optical Talbot 
carpet for monochromatic light to be scattered 
on 4-slit grating. A significant parameter in this 
optical pattern is the Talbot Length 



where d is the period of the diffraction grating 
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and A is the wavelength of the light incident 
on the grating. We will deal with this param- 
eter often enough at representations of interfer- 
ence patterns in near-field regions. Instead of 
monochromatic light scattering, however, we will 
simulate here particles (cold neutrons) scattering 
on nanoscale gratings (the period d is multiple 
of the particle wavelength). 
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FIG. 1: The optical Talbot effect for monochromatic light, 
shown as a "Talbot Carpet". The figure has been captured 
from http: / / en.wikipedia.org/wiki /Talbot_effect . 

On the other hand, as a distance from the slit 
grating to a detector increases the exotic near- 
field interference transforms to a diffraction pat- 
tern observed in the far-field region. Beautiful 
magnificence of the fractal structures disappears. 
Instead of it divergent rays from the slit source 
come into being. Intensity of these rays is de- 
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scribed by a well-known formula [10] 
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Function io(C) is an envelope describing diffrac- 
tion on a single slit and 
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is a phase shift of the waves emitted from two 
nearest slits, see Fig. |2} Observe that Eq. (J2| is 
quite common formula. For example, it describes 
also revolution of spin of neutrons flying through 
an iV-periodic magnetic field [TT] . 




FIG. 2: Diffraction from two slits in far-field region: d is 
distance between the slits, 9 is deflection angle with respect 
to axis z, and £ is phase shift of wave rays emitted from the 
slits. 

Most striking observation in Fig. 1 is existence 
of fractal structures, sizes of which are smaller 
than distance between slits. As was mentioned 
above, emergence of the fractal structures have 
been studied by Berry et al. [I], Amanatidis et 
al. [5], Sanz [6]. Sanz, in particular, has de- 
voted the studies whether one can depict by the 
Bohmian trajectories the quantum fractal struc- 
tures. The authors note that universal fractal 
features of quantum theory might be useful in 
the field of quantum information, for creating 
efficient quantum algorithms. 

As was mentioned, interference effects in the 
near-field regions can be important with the 
point of view of quantum computing perspective. 
Theoretical study of the near-field interference 
by means of preparation of quantum mechanical 
models is the first step for an understanding of 
the quantum computing perspective. Gaussian 



wavepacket [12], [TS] is a simplest model for study- 
ing iV-slit grating interference. Sec. [TT] deals with 
the Gaussian wavepacket and its Fourier trans- 
forms that give rise to emergence of complex, 

In Sec. 



Ill 



we sim- 



time-dependent, variance, 
ulate interference by radiation of the Gaussian 
wavepackets from the iV-slit grating. Their dis- 
persion in the near-field region produces an inter- 
ference pattern that manifests itself in a fractal 
organization of probability density of the wave 
function. As the detector is shifted in the far- 
field region the interference pattern transforms 
to diffraction pattern described by Eq. (|2]). In 
that region fractalality disappears. Instead of 
this, principal maxima of radiation come into 
being. And they are partitioned by (N — 2) sub- 
sidiary maxima. Next, the problem is to under- 
stand, how does a particle pass through the N- 
slit grating up to a detector. David Bohm had 
revealed that the particle travels along a single 
optimal trajectory [HIE], which is named in lit- 
erature the bohmian trajectory. This approach 
together with calculation of trajectories is dis- 



cussed in Sec. [TV] Emergence and development 
of the fractal Talbot patterns are considered in 



Sec. V at the slit grating containing many slits 
and at varying the period d of the grating. Com- 
parison of classical trajectories and bohmian is 



given in Sec. [VI] concluding section. In this sec- 
tion we discuss the bohmian trajectories passing 
through a slit up to the detector, and waviness 
of the bohmian trajectories caused by exchange 
of virtual particles with vacuum that is tuned by 
de Broglie pilot- wave. 




FIG. 3: Interference experiment in a cylindrical geometry - 
the grating is placed in the plane (x, y) and slits go along the 
axis y. Gray arrow points direction of wave radiation. 
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II. GAUSSIAN WAVEPACKET AND ITS 
FOURIER TRANSFORMS. 

As seen in Fig. |2j Eq. ^ results from com- 
putation of interference of cylindrical waves di- 
vergent from the slits. The waves in the vicin- 
ity of the slits are no cylindrical, in general. 
They have a complex form when transforming 
from a plane wave to cylindrical one. For the 
sake of simplicity, we suppose that slit borders 
are fuzzy under Gauss distribution, and a wave 
from the slit is approximated by the Gaussian 
wavepacket [12j [16] 
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Observe that v 9 o( x ) = P( x ) * s the probability 
density distribution that has a mean xq and vari- 
ance cr 2 > 0. 

We believe that in the immediate neighbor- 
hood on each slit, see Fig. |3j a wave field is 
described by the Gaussian wavepacket (|3]). It 
rely on assumption that edges of the slits are 
not ideal, but rather fuzzy. 

Observe that the wavepacket ^ disperses as 
it moves away from the slit to a region pointed 
out by grey arrow in Fig. |3j Let us express this 
wavepacket by superposition of harmonic waves 
exp{— 2Trikxo} with wave number k ranging from 
— oo to oo. To that end, we execute the Fourier 
transform 
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The function $o(&) is seen to be a harmonic wave 
exp{— 2nikxo} having an amplitude 

A{k) = V8^exp{-n 2 Acx 2 k 2 }. (5) 

As we move off from the slit, the wave com- 
ponents &o(k) spread along axis x differently 
for different k. Let the dispersed wave com- 
ponent be $>o(k) exp{— i27r&; ■ ttx} . Additional 
factor 7r scales the dispersion along axis x. Ob- 
serve that x = v t and a speed of the shift along 



axis x is v = Hk/m. Here H is the reduced 
Planck constant, m is a particle mass. As far 
as 2itk • k = (27r/c) 2 /27r, a shifted wave compo- 
nent is 

$ (A;)exp|-i(27rA;) 2 ^-t|. (6) 

Let us now execute the inverse Fourier trans- 
form of the function (|6l) 
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Here 



is a complex time-dependent spreading. It 

(8) 
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should be noted, that it is a complex variable 
of time t. Getting ahead, we can say that this 
complex variable [121 H3] determines fractal pat- 
tern in the near-field region. Observe that at 
t = we have cr t =o — °~ an d the function ^ 
comes to 4>q(x). 

III. RADIATION FROM iV-SLIT GRATING. 

The function Q is not yet a wave function 
since it contains no a term describing its trans- 
lation forward the region, as shown by grey ar- 
row in Fig. |3j The term describing such a trans- 
lation along axis z is represented by a factor 
exp{ico>t — ik z z}. Here E — tkv is a particle en- 
ergy and p z = hk z is its momentum along axis 
z. The wave function, in such reads 

xq, z) = cf)(x, xq, t) exp{iut — ik z z} (9) 
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A velocity along axis z is v z = p z /m = hk z /m. 
So, we can express time t in the function ^ 
via variable z/v z . Therefore, arguments in the 
function ^(x, x , z) do not contain t. 
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FIG. 4: Diffraction pattern in a transient region, (x;z) S 
((0 ■ ■ ■ 1.4 x 10" 6 ); (5 x 10" 8 ■ • ■ 3.8 x 10" B )) m, from near- 
field to far-field ones. Red arrows point out to directions of 
radiation along the principal peaks. Blue arrows point out to 
directions of weak radiation along the subsidiary maxima. 

Putting the frame of axis x in center of the 
slit grating, containing (N + 1) slits, (n = 



0, 1, 2, • • • , N), as shown in Fig. |3j we find that 
position of n th slit is Xq = (n — N/2)d. Here dis- 
tance d is period of the grating. Superposition 
of the waves ^ emitted by all (iV+1) slits reads 

\*{x,z)) = 7r!j^ 71 \ d < z 



N 



Probability density 



(10) 
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p(x,z) = (§?(x,z)\^(x,z)) 

is an observable. Fig. [4] shows the probability 
density p(x, z) which is calculated within a tran- 
sient region from the near-field region to the far- 
field ones. The grating contains four slits. It 
should be noted, that coordinate z in this figure 
begins from a value that is slightly more than 
zero. The whole point is that the wave function 
at z = has singularities located on the slits. 

Here we have simulated neutron scattering on 
the slits, with the neutrons having a wavelength 
A = 5nm. Kinetic energy of the neutrons is 
about 5.25 x 10~ 24 J. It means, that the neu- 
trons are very cold [33], i.e., temperature of the 
neutrons is about T = E/kB « 0.38 K. Here 
kB ps 1.38 x 10 -23 JK -1 is the Boltzmann con- 
stant. 
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FIG. 5: Diffraction in the far-field region, z = 0.004 m. 
Cross-section of the diffraction pattern (see previous figure) is 
displayed by circles. Solid curve shows the same cross-section 
calculated by the diffraction formula |2j| with parameters £ 
and 7o defined in Eqs. (l2| and (l3|. 



Fig. [4] demonstrates diffraction pattern in a 
transition region from the near-field region to 
the far-field one. Red arrows point out to the 
principal maxima. And blue arrows point out to 
the subsidiary maxima. Alternation of the prin- 
cipal maxima and subsidiary ones is described 
by Eq. Parameters ( and Jo, in this case, 
have the following forms pU] 
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where denominator D(z) is as follows 



D(z) = m 2 a 4 + h 2 (z/v z y/4:. 



(13) 



(14) 



Here (z/v z ) = t is a flight time along the path 
length z. Diffraction curve ( pj), at substituting 
terms ( and I by functions ( 12 )-( 13 ), has been 
calculated in the far-field region, z = 0.004 m, 
is shown in Fig. [5] Circles in this figure relate 
to results simulated by means of the Gaussian 
wavepacket. 
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FIG. 6: Diffraction pattern in the near-field region, (x, z) g 
((0 ■ • ■ 3 x 1CT 7 ); (8 x 10~ 10 • ■ • 2 x 10" 7 )) m. Ripples of the 
probability density lengthen as distance to the N-slit sources 
is increased. 

Turn back to Fig. |4} In order to show the 
diffraction pattern in the far-field region, we have 
coarsen resolution of the probability density dis- 
tribution. For this reason, we see a rough pile-up 
of maxima in the near-field region. Fig. [6] shows 
the same pattern at more detailed resolution in 
the near-field region. One can see, the proba- 
bility density distribution increases catastrophi- 
cally nearby the slits. At some distance from the 
slits, their radiations are superimposed with each 
other. Such a superposition manifests itself by 
emergent peaks. The peaks near the slits seem 
short. And they become longer, as a look moves 
from the slits to the far-field region. As soon as 



the detector is shifted to the far-field region, the 
peaks are transformed to typical diffraction rays 
going away to infinity. 

IV. DENSITY DISTRIBUTION PLOT AND 
BOHMIAN TRAJECTORIES. 

Let us project the density distribution p(x, z) 
to the plane (x,z). Fig. 7 demonstrates this 
projection in grey palette. Dark places show 
maximal values of p(x, z), black patches, in par- 
ticular, side with the slits. Light places relate 
to minimal values. Fractality in the interference 
pattern are well visible when passing from the 
slit sources to the right edge of the figure. 

Violet curves traced in the upper part of the 
figure depict bohmian trajectories. Finding the 
trajectories is based on variation of the action 
integral [XT] , which leads to the principle of 
least action. A general formula computing the 
bohmian trajectory [T31 [IB], the guidance equa- 
tion US [111, reads: 



x 



m 



*(x,z))- l V\*(x,z))) (15) 



As a result, we can find a current position of the 
particle by the following formulas: 



x(t) = Xq + v x dr 



(16) 



^ z(t) = z + v z ■ t. 



The velocity along z is v z = %k z /m. It stems 
from the term exp iut — ik z z as adopted before, 
see Eq. Q. 

Velocity v x is seen from Eq. (15) to be (a) 



proportional to gradient of the wave function; 
and (b) inversely proportional to the same wave 
function. It means: (a) a trajectory undergoes 
greatest variations in parts, where the wave func- 
tion demonstrates slopes; and (b) the trajectory 
avoids where the wave function tends to 

zero. Violet wavy curves in Fig. [7] manifest them- 
selves the above properties well enough. One 
can see, the wavy curves group predominantly 
in dark-grey spots and avoid white spots. In ex- 
treme cases, the trajectories traverse the white 
spots almost transversally. 
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FIG. 7: Density distribution plot drawn by grey palette displays the Talbot effect: white color relates to zero intensity, and 
black color points to maximal intensity. Violet curves drawn in the upper part of the plot depict the bohmian trajectories. 
Distance d = 5 x 10~ 8 m. It is equal to 10 • A, A = 5 x 1CT 9 m.. 



Does this bohmian trajectory pattern relate 
to the many-worlds theory? Answer is nega- 
tive [T9| [21] . For that aim, let us scatter particles 
on the N-slit grating by single-piece. Let a single 
particle pass, for example, the central slit nearby 
the top border. As is shown in Fig. [7j the particle 
scatters to a direction pointed out by blue arrow 
a. Next, the particle goes on a wavy stream and 
follows to a place pointed out by blue arrow b. 
Here we suppose that, movement of the particle 
submits to the principle of least action. There- 
fore, a number of problems arises right now. The 
first problem arising here is as follows: what is 
cause which forces the particle to change its own 
direction in the vicinity of point a? And the 



second problem is: what is cause which forces 
the particle to perform wavy motions? We con- 
firm that, any experimental setup for a quantum 
mechanical experiment is, in fact, the quantum 
instrument. The setup contains, in our case, N- 
slit grating with its parameters fitted for a wave- 
length of the particle. Boundary conditions of 
the setup determine the wave function on the 
edges. In whole, the experimental setup deter- 
mines configuration of the wave function given 
in the working space. The function is named de 
Broglie pilot-wave [T9J [20]. In fact, its squared 
image, the density distribution (11), is shown in 



Fig. [7] in grey palette. So, the bohmian trajec- 
tory shows an optimal path for a particle, which 
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FIG. 8: Interference from TV = 64 slits in the near-field region. Density distribution plot drawn by grey palette displays the 
Talbot effect: white color relates to zero intensity, and black color points to maximal intensity. Violet curves drawn in the 
upper part of the plot depict the bohmian trajectories. Grating parameters are A = 5 nm, d = 10A = 50 nm, zt = 1000 nm. 



is guided by the de Broglie pilot-wave within the 
working space. 



V. FRACTALS IN THE NEAR-FIELD REGION. 

In order to study fractal patterns arising in 
the near-field region in detail, amount of slits in 
the grating is not enough. In fact, the amount 
should tend to infinity. We will consider here, 
however, emission from a finite grating contain- 
ing 64 slits, see Fig. |8j One can see, in the vicin- 
ity of the slits placed in the middle of the grating 
(in the figure that place is drawn by red square) 
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FIG. 9: Diffraction in the far-field region, z — 1.25 m. Cross- 
section of the diffraction pattern is displayed by circles. Solid 
sharp principal peaks show the same cross-section calculated 
by the diffraction formula |2| with parameters £ and Iq de- 
fined in Eqs. (121 and Q 1 3 p . Dot envelope curve is Io{x) ■ N 2 . 



the Talbot carpets can be perfect enough. As 
the wave front spreads to an area reaching to 



the far-field region, the Talbot ordered struc- 
ture is dissolved by triangle-like manner, as is 
described in [12J. In the far-field region we will 
observe an usual diffraction from iV-slits, i.e., a 
set of principal maxima partitioned from each 
other by (N — 2) subsidiary maxima, Fig. [9j 

Violet wavy curves drawn in the upper part in 
Fig. [8] are Bohmian trajectories. They occupy 
preferably regions, where the probability den- 
sity reaches local maxima. And vice-versa, they 
avoid local minima. Instructive to compare be- 
havior of these Bohmian trajectories with those 
shown in Sanz's and Milet's article [12]. 
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FIG. 12: Talbot carpet: d = 40A = 200 run, z T = 16000 nm. 



A. Talbot carpets. 
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FIG. 10: Talbot carpet: d = 10A = 50 nm, z T = 1000 nm. 
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FIG. 13: Probability density distributions along cross-section 
pointed out by arrow a in Figs. |10| |11| and [12] are shown by 
green, red, and blue colors, respectively. Insert (a) shows 
behavior of p(x, z) near the origin of coordinates. 



Let us now discuss the Talbot patterns cut 
from the red square adjoining to the slit screen, 
see Fig. [8j Here we study emergence of the Tal- 
bot carpets at different relations of the period 
d and wavelength A. Since the slit grating con- 
tains only N = 64 slits, we can study a small set 
of such relations, d = kX, k is integer. We have 
given d = 10A, d = 20A, and d = 40A. According 
to Eq. (Jll), we have zt = 2kd. In the Talbot's 
units we can compare the Talbot carpets calcu- 



fig. 11: Talbot carpet: d = 20A = loo nm, z T = 4000 nm. lated for different input parameters. Figs. 10 11 



and 12 demonstrate the Talbot carpets calcu- 
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FIG. 14: Probability density distribution approaches infinite set of (^-functions as X/d tends to zero. 



lated at A = 5 nm and d = 10A = 50 nm, 
d = 20A = 100 nm, d = 40A = 200 nm, re- 
spectively. 

Instructive to compare Talbot carpets shown 
in Figs. [T] and 10 The both figures are seen to 
relate to each other as negative and positive im- 
ages. One can see, the both carpets show equiv- 
alent details. Whereas, the Talbot carpet shown 



in Fig. 11 is seen to contain more subtle details. 
In other words, this carpet is more fine-grained, 
than previous one. The Talbot carpet shown in 



Fig. 12 contains so much subtle details that the 
figure has been drawn only in the first quarter of 
the Talbot Length Zt- Here we can see fractal 
organization of the Talbot carpet very clearly. 

Behind the iV-slit grating, the field we can 



see in Figs. 10, 11 , and 12 has a fractal struc- 
ture of fantastic complexity. Let us glance on 
behavior of the probability density distribution 
along a cross-section pointed out by arrow a in 
these Figs. We see, that the probability den- 
sity vanishes at points z = and z = zt and it 
reaches maximal values at z = Zt/4, z = zt/2, 
z = 3zt/4. In the other points the probabil- 
ity density alternates minimal and maximal val- 
ues by irregular manner. More strictly, the ir- 



regularity exposes fractal nature. Fig. 13 shows 
behavior of the probability density distribution 
along the mentioned cross-section. Insert (a) 
shows its behavior near the origin of coordinates, 
z G (0,zt/10). One can guess, that the proba- 
bility density distribution tends to Cantor-like 
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set, as d/X tends to infinity. This is in good 
agreement with a statement given by Berry et 
al. in [3], that proclaims that the Talbot fractal 
emerges as the ratio X/d approaches zero and 
when the number N of illuminated slits tends 
to infinity as well [21 E]. At N finite, however, 
the Talbot patterns are blurring and defocusing. 
And what is more, the Talbot effect is destroyed 
catastrophically as an observer shifts the detec- 
tor either to edge of the iV-slit grating or beyond 
the near-field region. 

The Cantor set is a set of points lying on a sin- 
gle line segment that has a number of remark- 
able and deep properties [33]. First, the Can- 
tor set cannot contain any interval of non-zero 
length. On the other hand, integral of the prob- 
ability density distribution throughout all phys- 
ical space has to be equal to unit. From here 
it follows, that the probability density distribu- 
tion approaches an infinite set of 5-functions, see 



Fig. 14, as the ratio X/d tends to zero. 

Physically, limit of A — > is not available. It 
should be noted here, that at the wavelength 
tending to zero, kinetic energy of the incident 
particles approaches infinity. In that case, the 
particle beam will heat up the grating. And sec- 
ondly, if the kinetic energy is increased further, 
the particles begin to destroy the grating. 

In the upper parts of Figs. [TUJ [TTJ and [12 
Bohmian trajectories representing particle's 
paths have been drawn. They are pictured by 
violet dots tracking predominantly along dark 
places relating to heightened values of the prob- 
ability density. Their behavior is well visible in 
the vicinity of the dark nodes localized at cross- 
sections Zi/A, zt/2, 3zt/4. In general, particles 
jump along separated points until they leave the 
near-field region. Real possibility is that, the 
particle is tunneling throughout the suppressed 
intervals. In other words, it reaches the far-field 
region along fantastical zigzag paths, that are 
frequently interrupted by tunneling. One can 
see, behavior of the Bohmian trajectories is com- 
plex enough in the fractal media. Nevertheless, 
it is predictable, since it is based on solution of 
two coupled equations - the quantum Hamilton- 
Jacobi equation and the continuity equation. 
The both result from the Schrodinger equation. 



VI. CONCLUDING DISCUSSION. 

A fractal, as defined by Mandelbrot, "is a 
shape made of parts similar to the whole in 
some way" [22] • An exact fractal is an "object 
which appears self-similar under varying degrees 
of magnification in effect, possessing sym- 
metry across scale, with each small part replicat- 
ing the structure of the whole" [22J. In this key 
Berry has also written in a private communica- 
tion the following assertion "Fractality requires 
three conditions: infinitely many slits, paraxial 
propagation, and discontinuous initial conditions 
(i.e. sharp slits). Then the fractal structure, in- 
cluding fractal dimensions, can be calculated ex- 
plicitly." As for the fractal Talbot effect, Berry 
et al. have written in [1] "Quantum and optical 
carpets provide a dramatic illustration of how 
limits in physics that seem familiar can in fact 
be complicated and subtle. It is no exaggeration 
to say that perfect Talbot images, and infinite 
detail in the Talbot fractals, are emergent phe- 
nomena: they emerge in the paraxial limit as 
X/d (here d is signed instead of a, V. S.) ap- 
proaches zero. At first this seems paradoxical, 
because the short-wave approximation is usually 
regarded as one in which interference can be ne- 
glected, whereas the Talbot effect depends en- 
tirely on interference. The paradox is dissolved 
by noting that the Talbot distance increases as 
the wavelength approaches zero, so here we are 
dealing with the combined limit of short wave- 
length and long propagation distance: in the 
short-wavelength limit, the Talbot reconstructed 
images recede to infinity." 

Talbot patterns arising in the near-field region 
demonstrate the above mentioned signs fully. 
These patterns have emergent at simulation of 
scattering cold neutrons on many-slit grating (we 
have chosen parameters relating to this particle). 
A question arises here - how does a single parti- 
cle travel behind the grating? 

First of all, let us recall some solutions from 
the classical mechanics [23J. 

Geodesic trajectories in classical mechanics to 
be found by applying the principle of least ac- 
tion [23] are real trajectories of mechanical ob- 
jects. A beam of the trajectories, that tracks a 
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relief, stems from the initial conditions slightly 
differing from each other. Observe that the 
geodesic trajectories nowhere never intersect. 

The number of trajectories, crossing a surface 
8S, is conserved, no matter how the surface de- 
forms at moving along the beam. In fact, this 
observation demonstrates the conservation law 
of number of the trajectories passing through 
the surface, i.e., the trajectories never disappear 
and does not appear again. Violation of this law 
in classical physics would be tantamount to rec- 
ognize teleportation - a classic body disappear 
suddenly and unexpectedly announced again. 

Interpretation of the continuity equation of the 
density of trajectories is guided by incompress- 
ible fluid, which flows along routes specified by 
the geodesic trajectories. All physical space we 
can imagine is filled with this fluid [23j 121] . The 
basis for such an idealization is experience which 
shows, that there is a rather broad class of fluids 
for which even large changes in pressure do not 
lead to significant change in density. This fluid 
fills the environment continuously. Its molecu- 
lar structure, at that, is ignored. In classical 
physics the continuity equation does not deter- 
mine the fate of trajectories in the future. Their 
fate is determined only by the Hamilton- Jacob i 
equation. Other situation arises in quantum me- 
chanics. Here the both equations, the conti- 
nuity equation and the Hamilton- Jacobi equa- 
tion, connected with each other via the bohmian 
quantum potential [HI EE], take part in determi- 
nation of the geodesic trajectories. First scientist 
was Madelung who had derived the same set of 
equations in 1926 [25] . 

In accordance with the de Broglie-Bohm inter- 
pretation, the wave function sets up a ^-field, 
which fills physical space of the experimental 
setup. Its squared modulus, the probability den- 
sity distribution, is shown in Figs. [7|8p0}|l2t in 
gray palette. All physical space we can guess is 
filled by a fluid-like background with the density 
distribution |\1>| 2 [25, 26J. It is important to note, 
that the state of such a fluid is determined by ge- 
ometry of the experimental setup. And route of 
the trajectories depends on the density of the 
fluid, its gradients, in the neighborhood of each 



point of the experimental setup. The guidance 
equation allows to find the trajectories that pen- 
etrate the field by the best way. In turn, the 
density distribution depends on the route of the 
trajectories. They are mutually dependent. 

Such a fluid is seen to be vacuum. Polarization 
of vacuum is an extraordinary phenomenon since 
Hendrik Kasimir had proposed the existence of 
physical forces arising from a quantized field • 
From this perspective, change of direction of a 
trajectory of a particle can be expressed in terms 
of exchange of virtual particles in vacuum. The 
vacuum has a vastly rich structure. It, implicitly, 
has all of the properties that a particle may have. 
This situation can be expressed by a sentence - 
the vacuum contains relative ' 'nothing 1 , and at 
the same time, the potential ' all 1 . 

The Feynman path integral [2E] is a way to 
understand many manifestations of the vacuum. 
In fact, a proposition is as follows: virtual pairs 
emergent in a short time and annihilating in the 
end of the time loop are pairs permitting the 
particle to test a trajectory both forward and 
backward in time [291 USB EI]- The extra de- 
gree of freedom obtained by allowing both direc- 
tions in time gives the particle to get a non-local 
information needed for establishing an optimal 
path from a source to detector without the in- 
volvement of intelligent particles, intelligent ob- 
servers, or multiple universes. 



Acknowledgments 

The author thanks Dr. A. S. Sanz, who im- 
pelled more deep study of fractal interference 
patterns in the near-field region. The author 
thanks also Prof. M. V. Berry for valuable com- 
ments relating to the fractal Talbot effect. The 
author expresses his sincere thanks to Miss Pipa, 
the site administrator of Russian Quantum Por- 
tal, for developing and writing a program calcu- 
lating the density distribution of the wave func- 
tion at scattering particles on iV-slit grating. 
The program draws also Bohmian trajectories 
outgoing from the slits. 



12 



A. D. Cronin, J. Schmiedmayer and D. E. Pritchard. 
Atom Interferometers, e-print URL http : // arxiv . org/ 
abs/0712 3703vl| (21 Dec 2007). 

M. V. Berry and S. Klein. Integer, fractional and frac- 
tal Talbot effects, JOURNAL OF MODERN OPTICS 
43(10) (1996) 2139-2164. 

M. V. Berry. Quantum fractals in boxes, J. Phys. A: 
Math. Gen. 29 (1996) 66176629. 

M. Berry, I. Marzoli, and W. Schleich. Quantum carpets, 
carpets of light, Phys. World (6) (2001) 1-6. 
E. J. Amanatidis, D.E. Katsanos and S.N. Evan- 
gelou. Fractal Noise in Quantum Ballistic and Diffusive 
Lattice Systems, e-print URL http://arxiv.org/abs/ 



quant-ph/0310125 ( (20 Oct 2003). 
A. S. Sanz. A Bohmian approach to quantum fractals, J. 
Phys. A: Math. Gen. 38 (2005) 60376049. 
J. F. Clauser and J. P. Dowling. Factoring integers with 
Young's N-slit interferometer, PHYSICAL REVIEW A 
53(6) (1996) 4587-4590. 

A. Steane. Quantum Computing, Rept. Prog. Phys. 61 
(1998) 117-173. 

S. Olmschenk, D. N. Matsukevich, P. Maunz, D. Hayes, 
L.-M. Duan, and C. Monroe. Quantum Teleportation 
Between Distant Matter Qubits, e-print URL http:// 



arxiv . org/abs/0907 . 5240 (30 Jul 2009) 



V. I. Sbitnev. Particle scattering on a N-slit screen: Tal- 
bot carpet and far-field diffraction, Kvantovaja magija 
6(1) (2009) 1101-1112, in Russian. 

M. M. Agamalyan, G. M. Drabkin, and V. I. Sbitnev. 
Spatial spin resonance of polarized neutrons. A tunable 
slow neutron filter, Phys. Rep. 168(5) (1988) 265-303. 
A. S. Sanz and S. Miret-Artes. A causal look into 
the quantum Talbot effect, J. Chem. Phys. 126 (2007) 
234106; e-print URL http://arxiv.org/abs/quant-ph/ 
0702224, (19 Jun 2007). 

A. S. Sanz and S. Miret-Artes. A trajectory-based un- 
derstanding of quantum interference, J. Phys. A: Math. 



Gen. 41 (2008) 435303; e-print URL http://arxiv.org/ 



I abs/08 06.2105 (1 Oct 2008) 
[14] D. Bohm, A suggested interpretation of the quantum the- 
ory in terms of "Hidden Variables", Phys. Rev. 85 (1952) 
166-193. 

D. Bohm and B. J. Hiley. The Undivided Universe: An 
ontological interpretation of quantum theory, (Routledge, 
London, 1993). 

G. I. Mark. Analysis of the spreading Gaussian 
wavepacket, Eur. J. Phys. 18 (1997) 247-250. 
V. I. Sbitnev. Bohmian trajectories and the path integral 
paradigm. Complexified Lagrangian mechanics IJBC, 
19(7) (2009) xxxx-x xxx: e-print URL [http: //arxiv .| 
I org/abs/0808.12"45l (8 Aug 2008). 

[18] M. Davidovic, D. Arsenovic, M. Bozic, A. S. Sanz, and S. 
Miret-Artes. Should particle trajectories comply with the 
transverse momentum distribution?, Eur. Phys. J. Spe- 
cial Topics 160 (2008) 95, e-print URL [http: //arxiv .| 
I org/abs/0803 .26061 ( 18 Mar 2008). 

[19] A. Valentini. De Broglie-Bohm Pilot- Wave Theory: 
Many Worlds in Denial? In: Everett and his Critics, eds. 
S. W. Saunders et al.. (Oxford, University Press, 2009), 
e-print URL |http : //arxiv ."or g/abs/08 11 . 0810"] (5 Nov 
2008). 



[20] W. Struyve and A. Valentini. de BroglieBohm guidance 
equations for arbitrary Hamiltonians. J. Phys. A: Math. 
Theor. 42 (2009) 035301 (18pp). 

[21] A. Kent. One world versus many: the inadequacy of 
Everettian accounts of evolution, probability, and sci- 
entific confirmation, In Many Worlds? Everett, Quan- 
tum Theory and Reality, eds. S. Saunders, J. Barrett, 
A. Kent and D. Wallace, O.U.P. upcoming, e-print URL 
|http : //arxiv . org/abs/0905 ■ 0624| (5 May 2009). 

[22] P. S. Addison, Fractals and Chaos (Institute of Physics, 
Bristol, 1997). 

[23] C. Lanczos. The variational principles of mechanics, 

(Dover Publ., Inc., N. Y. 1970). 
[24] R. E. Wyatt. Wyatt2005 Quantum Dynamics with Tra- 
jectories: Introduction to Quantum Hydrodynamics, 
(Springer, Berlin, 2005). 
[25] E. Madelung. Quantentheorie in Hydrodynamischer 

form. Zts. f. Phys, 40 (1926) 322 326. 
[26] D. Bohm and J. P. Vigier. Model of the Causal Inter- 
pretation of Quantum Theory in Terms of a Fluid with 
Irregular Fluctuations, Phys. Rev. 96(1) (1954) 208 - 216. 
[27] H. G. B. Casimir. On the attraction between two per- 
fectly conducting plates, Proc. Con. Ned. Akad. van 
Wetensch B51(7) (1948) 793-796. 
[28] R. P. Feynman and A. Hibbs. Quantum Mechanics and 

Path Integrals, (McGraw Hill, N. Y. 1965). 
[29] G. N. Ord and R. B. Mann. Entwined pairs and 
Schrodinger's equation, Annals of Physics, 308(2003) 
478-492; e-print URL [http : //arxiv . org/abs/ 



quant-ph/0206095, (14 Jun 2002). 
[30] G. N. Ord, J. A. Gualtieri, and R. B. Mann. A physi- 
cal basis for the phase in Feynman path integration, e- 
print URL |http : //arxiv . or g/abs/quant-ph/041 1005 
(31 Oct 2004). 

[31] P. J. Werbos and L. Dolmatova. The Backwards- 
Time Interpretation of Quantum Mechanics - Revisited 
With Experiment, e-print URL http://arXiv.org/abs/ 
quant- ph/0008036, (7 Aug 2000). 

[32] http://en.wikipedia.org/wiki/Talbot_effect 

[33] http:/ /en. wikipedia.org/wiki/Thermal_neutron 

[34] http://en.wikipedia.org/wiki/Cantor_set 



